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ABSTRACT 



Aims. The properties of the Galactic supernova remnant (SNR) RX J0852. 0-4622 are theoretically analysed. 

Methods. An explicitly time-dependent, nonlinear kinetic model of cosmic ray (CR) acceleration in SNRs is used to describe the 
properties of SNR RX J0852. 0-4622, the accelerated CRs and the nonthermal emission. The source is assumed to be at a distance of 
~ 1 kpc in the wind bubble of a massive progenitor star. An estimate of the thermal X-ray flux in such a configuration is given. 
Results. We find that the overall synchrotron spectrum of RX J0852. 0-4622 as well as the filamentary structures in hard X-rays lead 
to an amplified magnetic field B > 100/iG in the SNR interior. This implies that the leptonic very high energy (VHE) 7-ray emission 
is suppressed, and that the VHE 7-rays are hadronically dominated. The energy spectrum of protons produced over the life-time of the 
remnant until now may well reach "knee" energies. The derived 7-ray morphology is consistent with the H.E.S.S. measurements. The 
amount of energy in energetic particles corresponds to about 35% of the hydrodynamic explosion energy. A remaining uncertainty 
concerns the thermal X-ray flux at 1 keV. A rough estimate, possibly not quite appropriate for the assumed wind bubble configuration, 
results in it being larger than the nonthermal flux at this energy. 

Conclusions. It is concluded that this SNR expanding into the wind bubble of a massive star in a dense gas environment can be a 
hadronic 7-ray source that is consistent with all existing multi-wavelength constraints, except possibly the thermal X-ray emission. 
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1. Introduction 



-(-acceleration of particles 
- gamma-rays:theory 



shock waves - supernovae individual (SNR RX J0852. 0-4622) 



RX J0852.0-4622 (also known as G266.2-1.9) is a shell-type su- 
pernova remnant (SNR) with a diameter of 2°, located in the 
Galactic plane. The SN R was originally discovered in X-rays 
• £h ' with t he ROSAT satellite (lAschenbachlll998l;lAschenbach et al.L 
'. I1999I) . In projection along the line of sight, RX J0852.0-4622 
5_i ■ lies entirely within the still much larger Vela SNR and is only 
C$ \ visible in hard X-rays, where the thermal radiation from the 
Vela SNR is no longer dominant. While nonthermal emission 
from the shell of RX J08 52.0-4622 has been confirmed by sev- 
eral X-ray observa tories (ISlane et al.1 l2001b iBamba et all 120051; 
Ivud in et aTll2005l) . a clear detection of the thermal X-ray emis- 
sion from the shell or the interior was not yet possible because 
of confusion with the Vela SNR. This implies that the thermal 
emission is very wealfl 

The radio emission of RX J0852. 0-46 22 is also weak. In 
fact the SNR radio shell was only identified dCombi et all 1 19991; 
iDuncan & Greenl [2000; Stu par et all 120051) after its discovery 
in X-rays. Before that only a bright radio spot was known as 
"Vela Z" (Milne, 1968), which was usually identified with the 



Send offprint requests to: H.J. Volk 

* now at Institut fiir Astronomie und Astrophysik, Universitat 
Tubingen, Sand 1, 72076 Tubingen, Germany 

1 Very recently Uchiyama (2008) has indicated that the X-ray emis- 
sion from RX J0852.0-4622 might show a thermal component. 



Vela SNR. In addition, significant Galactic background varia- 
tion over the size of the remnant cannot be excluded. The radio 
spectrum of RX J0852. 0-4622 is therefore not well determined. 
Only for the northeaste rn rim a spectral index can be derived 
with moderate accuracy dDuncan & Greenl 12000). 

The shell of RX J0852.0-4622 was also detected in 
very high energy (VHE) 7- r ays by the H.E.S.S. collabora- 
tion (Ah aronian et all I2005L l2007al) . with a 7-ray flux at 
1 TeV as large as that from the Crab Nebula. Emission 
from the northwestern rim had be en detected al r eady before 
by the CANGAROO experiment dKatagiri et all 120051). The 
CAN GAROO data have been revised since then (Enomoto et all 
120061) . 

RX J0852.0-4622 i s the se cond SNR after RX J1713.7-3946 
(e.g. lAharonian et all H007b), where morphologically a SNR 
shell was unambiguously identified to accelerate particles to TeV 
energies and beyond. As a third and fourth source of this char- 
acter the objects RCW86 dHoppe & Lemoine-Goumardl 1 2008; 
l Aharonian et all 120081) and SN 1006 dNaumann-Godo et all 
2009) have recently been detected. A few other TeV sources 
have been detected in the H.E.S.S. Galactic plane survey that are 
spatially coincident with radio shell-type SNRs. So far, however, 
the data do not permit to unambiguously identify the 7-ray emis- 
sion with the SNR shells, e.g. using morphological arguments. 
In several cases, the TeV emission might also be associated with 
X-ray emitting pulsar wind nebula candidates. 
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Despite the scarcity of precise information from radio data 
and despite uncertainties about key astrophysical parameters, the 
prominence of RX J0852. 0-4622 has therefore prompted us to 
model the acceleration of both electrons and protons together 
with their nonthermal emission in detail, applying explicitly 
time-dependent nonlinear kinetic theory. The theory couples par- 
ticle acceleration on a kinetic level with the gas dynamical evo- 
lution of th e system dBerezhko et all ll996llBerezhko & Volkl 
119971 120001) . 

Compared to the other SNRs that were succe ssfully de- 
scribed within the framework of this theory (e.g. iBerezhkol 
120051 I2008t iBerezhko & Volkl 120061) . the present uncertain- 
ties regarding RX J0852. 0-4622 are quite large. Such important 
astronomical parameters as the distance and age are not well 
known. It is in fact not even clear, whether the source is in front 
or behind the Vela SNR. The latter object i s generally consi d- 
ered to he at a distance d = 250 ± 30 pc dCha et all [19990. 
A position behind t he Vela SNR cou ld correspond to a solution 
with d = 1 kpc dSlane et al. , 2001), whereas another solution 
could correspond to the earlier distance estimate of d = 200 pc 
dAschenbachlll998l) . 

This prompted us originally to construct indeed two quite 
different source scenarios, one in front, and the other behind 
the Vela SNR. They were to correspond to earlie r distance es- 
timat es: a "nearby " solution with d = 200 pc (lAschenbachl 
1998 ) in front of, and a "distant" solution with d = 1 kpc 
dSlane et all 1200 ll) behind the Vela SNR. However, no "nearby" 
solution could be found that fulfilled all the observational con- 
straints. Therefore we have abandoned the possibility of a small 
distance to RX J0852. 0-4622 and will assume for the sequel that 
the source is at a distan ce d = 1 kp c. 

In agreement with ISlane et ail (1200 lb but with independent 
additional arguments we shall assume that the observed non- 
thermal emission of RX J0852.0-4622 indicates that the SNR 
emerged from a core collapse explosion into the wind bubble of 
a massive progenitor star in a dense gas environment, possibly 
a molecular cloud. In this case the major part of the swept-up 
volume has originally been occupied by the highly diluted bub- 
ble gas that also has a minimal thermal emissivity. At the current 
epoch, however, we assume that the SNR shock already propa- 
gates into the dense shell of ambient interstellar medium (ISM) 
which has originally been compressed by the stellar wind. This 
also implies that the magnetic field upstream of the SNR shock 
is of interstellar origin. 

We note here that this solution has a similar character as 
the solution sugge s ted ea rlier for the object RX J17 13.7-3946 
dBerezhko & Volkl 120061) . Indeed, in this sense the two SNRs 
RX J0852.0-4622 and RX J1713.7-3946 can be considered as 
twins. 

The hydrodynamic state of the system, i.e. the given linear 
radius for known angular radius and distance together with the 
present age and shock velocity, etc., is basically determined by 
the choice of mechanical explosion energy, ejected mass, and 
external density. 

The lack of knowledge of the radio spectral shape makes it 
impossible to derive from the radio synchrotron observations - 
and as far as the magnetic field is concerned, from the combined 
radio and X-ray synchrotron spectral observations - the most 
relevant pair of physical parameters for the acceleration theory, 
namely th e proton injection rate and th e effec t ive magnet i c field 
strength dBerezhko et all 120021 IVolkl l2004t IBerezhkol 120051 

2 The distance to the Vela pulsar is estimated at d = 2871 ^pc 
dDodson et alll2003h . see also lCaraveo et alll l200lh . 



2008). Thus, even if we consider the circumstellar medium 
(CSM) structure to be given, it is not possible to directly pre- 
dict the form of the overall synchrotron and the full VHE 7-ray 
emission from theory. 

The observed overall nonthermal spectral shape - including 
the VHE 7-rays - and the small-scale filamentary structures in 
the nonthermal X-ray emission of RX J0852. 0-4622 neverthe- 
less provide evidence for effective nuclear cosmic ray (CR) ac- 
celeration, associated with considerable magnetic field amplifi- 
cation. This conclusion is possible because the overall nonther- 
mal spectrum can be theoretically fitted with an appropriate pro- 
ton injection rate, electron-to-proton ratio, and effective mag- 
netic fie ld strength (assum e d to be uniform inside the shocked 
CSM cf. IBerezhko & Volkl d2004b1) ). 

The main result of this paper is that the resulting VHE 7-ray 
flux is hadronically dominated. The well-known difference in the 
effectiveness of the basic radiation mechanisms then implies that 
the energy density of the nuclear component of the nonthermal 
charged-particle population in the SNR by far dominates that of 
the energetic electron component generated in the source. The 
energy in nonthermal particles at the present epoch amounts to 
~ 35 percent of the assumed total mechanical energy E sn — 
1.3 x 10 51 erg, released in the explosion. Therefore, from the 
point of view of energetics, this solution for RX J0852. 0-4622 
more than fulfils the average requirement on a SNR source of the 
Galactic CRs. In this context we shall also qualitatively discuss 
the question of possible escape of the highest-energy particles, 
accelerated at an earlier phase of the SNR than the one we can 
observe at present. 

In the next section the theoretical model is described. Section 
3 presents our assumptions regarding the values of the physical 
parameters as they are suggested by the broadband data and by 
physical considerations. It contains also a discussion of the ther- 
mal emission, even though it has not been possible to estimate it 
appropriately for the assumed wind bubble and shell configura- 
tion. The results for the gas dynamics, the particle acceleration 
coupling with it, the 7-ray emission, and for the thermal X-ray 
emission are presented and discussed in section 4. In section 5 
our conclusions are summarized. 



2. Model 

The theoretical model for the particle acceleration combined 
with the gas dynamics of the explosion has been described ear- 
lier, fo£instenc^Jn_are£ent analysis of SNR RX J17 13.7-3946 
(see IBerezhko & V5ikT l2006. and references therein). 

The ejected mass M c] has initially in its fastest moving outer 
parts a power law distrib ution dM^/du o c u 2 ~ k in flow vel ocity 
u, with 7 < k < 12 dJones et all fEtlt IChevalierl fl982h . We 
shall choose k = 8 (as in the case of SNR RX J1713.7-3946). 

The interaction of the ejecta with the CSM creates a strong 
outward-propagating shock wave in the CSM at which particles 
are accelerated. Our nonlinear model is based on an explicitly 
time-dependent solution of the CR transport equations together 
with the gas dynamic equations in spherical symmetry. In par- 
ticular the theory takes into account the adiabatic energy losses 
of thermal and nonthermal particles in the SNR interior, the dif- 
fusion of nonthermal particles from the outer shock into that in- 
terior, and the backreaction of CRs on the shock structure and 
dynamics. This backreaction decelerates the thermal gas already 
in front of the shock and leads to a smooth shock precursor that 
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reduces the Mach number of the subsequent collisionless plasma 
shock (the subshock), heating the inflowing gafl 

However, nuclear particles are only effectively injected into 
the acceleration process at those parts of the moving shock sur- 
face, where the locally mean magnetic field vector is quasi- 
parallel to the shock normal. These regions are characterized by 
magnetic flux tubes within which the injection of moderately 
suprathermal particles proceeds. 

In quasi-perpendicular shock regions, on the other hand, in- 
jection is instantaneously strongly suppressed or completely pro- 
hibited. In the extreme there are then separated regions at the 
shock surface, where particle injection is permanently either al- 
lowed or prohibited. Then also no nuclear particles are accel- 
erated in the prohibited regions. An example for this extreme 
situation is SN 1006, where the X-ray emitting polar cap re- 
gions correspond to the quasi-parallel regions. In these allowed 
regions diffusive shock acceleration strongly proceeds on the 
Bohm diffusion level (see below) as a result of the effective 
injection of low-energy particles. However, since Bohm diffu- 
sion is isotropic, the energetic particles can also cross field lines 
and thus possibly reach quasi-perpendicular regions of the shock 
where they can also accelerate. The extent to which this happens 
depends on the spatial scales that separate quasi-perpendicular 
regions from quasi-parallel regions. For a homogeneous exter- 
nal magnetic field about 80% of the shock surface is quasi- 
perpendicular in the above sense, and therefore onl y a fraction 
/re ~ 0.2 of the shock is efficiently accelerating (IVolk et ail, 
2003). However we shall argue in section 3.1.4 that for a SNR, 
propagating into a stellar wind bubble with a radiatively cooling 
shell of high-density gas, the above spatial scales are probably so 
small that / re w 1. This has substantial implications especially 
for the resulting thermal X-ray emission, because the observed 
VHE 7-ray emission (which is by implication of hadronic ori- 
gin) then requires a lower gas density. 

For given magnetic field strength the electron injection rate 
in spherical symmetry is determined by the intensity of the ob- 
served overall synchrotron spectrum. For given injection rate of 
nuclear particles and magnetic field strength the ratio K ep be- 
tween the densities of nonthermal electrons and nuclear particle 
can then be calculated. 

As a result of the streaming instability the accelerating CRs 
very effectively excite l arge-amplitude magnetic fluctuations up- 
stream of the SN shock (l^fT^lBlandford & OstrikerlfT97l 
iMcKenzie & Volkl 119821) . Since these fluctuations scatter CRs 
extremely strongly, we approximate the CR diffusion coefficient 
k(j>) by its lower limit, corresponding to a scattering mean free 
path equal to the particle gyro radius. In this so-called Bohm 
limit k(p) = [mc 3 /(3eB)](v/c)(p/mc), where e and m are 
the particle charge and mass, v and p denote the particle ve- 
locity and momentum, respectively, B is the effective magnetic 
field strength (see below), and c is the speed of light. Regarding 
the nuclear particles with the highest energies this Bohm limit 
may imply an underestimate for n(p) also for another reason, 
since for these particles the effective, amplified fiel d has spatial 
scales that are smalle r or equal to their gyro radius dBelll [2004; 
iPelletier el al.L 120061) . To this extent our assumption of Bohm 
diffusion in the amplified field yields an upper limit to the maxi- 
mum energy of the nuclear particle (e.g. Ziraka shvili & Ptu skin. 
2008). In addition we assume that the interior effective magnetic 



3 It is implicitely assumed here that any energy loss due to radiative 
gas cooling has no effect on the shock structure. Given the rather high 
internal temperature T su b of the remnant at this stage (see section 4.2) 
this appears to be a safe assumption. 



field has a uniform strength after its MHD compression in the 
thermal subshock. Practically speaking we assume this unifor- 
mity over a spatial scale that is la rge compared to the thick ness 
of the observed X-ray filaments (Bere zhko & VolkL l2004b). For 
a different point of view, see lPohl et al.l d2005l7 ! 

Another important nonlinear effect of the strong excitation 
of magnetic fluctuations by the accelerating particles themselves 
is the heating of the thermal plasma in the shock precursor that 
is generated by the accelerating particles. Combined analyti- 
cal and numerical efforts, using plasma theory t o give a non- 
linear d escription of the magn e tic field evolution [Luc ek & Bell 
(120001): [Bell & Lucekl 12000); iPtuskin & Zirakashvilil d2003l): 
iBelll d2004 : IPelletier el all d2006l) . concluded that a considerable 
amplification of the upstream magnetic field should occur in the 
acceleration process to what we call the effective magnetic field. 
The physical reason is that the beam of efficiently accelerated 
nuclear CR component excites in particular also a non-resonant 
magne tic mode in addition to the well-known resonant Alfven 
waves del [2004). The latter have nevertheless been argued to 
contribute dominantly to the o verall turbulent magnetic energy 
density in the shock precursor dPelletier el all 120 06). However, 
the three-di mensional M HD simulations of the non-r esonant in- 
stability bv lBelll d2004l) and lZirakashvili et all d2008l) show that 
the nonlinear growth of the magnetic fluctuations is accompa- 
nied by the formation of internal shocks and correspondingly 
strong dissipation which heats the thermal plasma^. Analogous 
dissipation should occur in three dimensions for the wave modes 
of the resonant streaming instability. We approximate this phys- 
ical process by assuming that the heat input into the thermal gas 
equals the (linear) growth of the turbulent field energy in the 
excited Alfven waves in the alread y amplified effective field (see 
[Berezhko, 20081: 1 Volk et aU l2008. for the correspondence of this 
approach to existing theory and experiment). 

As already mentioned in the Introduction, without a reliable 
spectral index for the observed spatially-integrated radio syn- 
chrotron emission we shall choose the strength of the effective 
field as well as the nuclear injection rate 77, such as to opti- 
mally fit the calculated synchrotron spectrum (from radio to X- 
ray energies) together with the calculated VHE 7-ray spectrum 
to the observations. Subsequently this spectrally fitted magnetic 
field is compared with the field derived from the observed fila- 
mentary structure. The degree of agreement between these two 
field strengths is then used as a measure of the success and self- 
consistency of the model. 

The filament-based magnetic field strength B& downstream 
of the shock is determined by the observed width L of the X-ray 
filament - interpreted as the synchrotron cooling length behind 
the shock - through the relation 

B d = [3 m y/(4erlll)}^(VTTP-S)-y, (1) 

where 5 2 = 0.12c/ '(r v)[V s /{ac)} 2 , l 2 « L/7 is the radial 
width of the X-ray emissivity q v (e u , r), e v is the X-ray energy, 
corresponding t o the frequency v, and rn denotes the classical 
electron radius (Berezhko & V61k, l2004al) . This field strength is 
clearly a lower limit. 

Since we have already used the observed amplitude of the 
VHE 7-ray spectrum to determine 77 we do not predict this am- 

4 Particle-in-cell simulations by iNiemiec et alj d2008l) led these au- 
thors even to the extreme conclusion that the non-resonant instability 
amplitudes become never large as a result of the bulk acceleration of 
the thermal gas by the streaming CRs. We believe that in a quasi-steady 
shock configuration the CR current is being steadily driven through the 
upstream gas by the downstream overpressure, as implied above. 
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plitude, as one could do using a detailed knowledge of the inte- 
grated synchrotron spectrum. The above consistency condition 
for the effective field still needs to be fulfilled for the solution to 
be acceptable. 

Moreover, we shall not only investigate whether the chosen 

values for B and ri are consistent with the Chandra measurement 

r 1 1 1 

of the X-ray filamentary structure (Ba mba et all 2005), but also 
to wh ich extent they are consi stent with the semi-empirical rela- 
tion dBerezhko & Valid 120061) 



Table 1. Key model parameters, and corresponding spectral, dy- 
namical, and morphological values expected from the calcula- 
tions. 



B$/(8nP c ) w 5 x 10" 



(2) 



which connects the upstream magnetic field pressure Bq / (8tt) in 
the shock precursor (i.e. upstream, but amplified by the CR insta- 
bility) with the pressure P c of the accelerated particles that drives 
the field amplification in the first place. Eq.(2) holds for a num- 
ber of SNRs that could be analyze d with the aid of a well-known 
radio spectrum for these sources (IVolk et al.Ll2005h . The de gree 
of fulfilment of this relation is a further quality criterium for the 
model to judge its success in theoretically describing the TeV 7- 
ray so urce. Eq.(2), or a relation of a similar type ( Be ll & LucekL 
2001), is likely to hold for an individual object also during its 
time evolution. Yet, in order to avoid the introduction of a fur- 
ther theoretical parameter, we shall consider Bo as constant in 
time, equal to its present value, in our evaluation of the models 
for RX J0852. 0-4622. We shall come back to this point in section 
4.3. 

In this specific form the three theory "parameters" B, r\, 
and K ep are determined quantitatively by comparison with 
the observations at the present age of the source. Their time- 
dependence during the evolution of the SNR is disregarded in 
this paper. 

The numerical solution of the dynamical equations at each 
instant of time yields the CR spectrum and the spatial distribu- 
tions of CRs and thermal gas. This allows the calculation of the 
spectra of the expected fluxes of nonthermal emission produced 
by the accelerated CRs, the morphology of the emission, and the 
future evolution inasmuch the same physical processes continue 
to work at later times. 

In the following we shall consider the wind bubble scenario 
for RX J0852. 0-4622 in the general framework of this model. 



3. Physical parameters of RX J0852.0-4622 

In this section, we will describe the physical parameters of 
the model. Section 3.1 gives an overview over the available 
broadband data. Section 3.2 describes the setup of the scenario. 
Because of the scarcity - and sometimes ambiguity - of the 
available data, we will use some general arguments from the 
non-thermal particle emission not only to constrain the accelera- 
tion parameters but also the environmental parameter gas density 
and its spatial distribution. The values of all relevant physical pa- 
rameters are given in Table 1 . 



3. 1. Broadband data of RX J0852.0-4622 

3.1 .1 . Morphology of RX J0852.0-4622 in general 

RX J0852.0-4622 is a shell-type SNR with a shell diam- 
eter of 2°, as seen in hard X-rays jAschenbachl Il998t 
lAschenbach et all 1 1 999h ISlane et all |2001|) and VHE 7-rays 
dAharonian et alir2005ll2007a) . Also the radio continuum emis- 
sion correlates spatially well with th e high energy emission 
dDuncan & Greenl 120001: IStupar et all 120051; lAharonian etaE 
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Bd from X-ray filaments 


139 /zG 


-Esn 

N S (R S ) 
7V g (r = 0) 
k 

Bo 
V 

K cp 

frc 


1.3 x 10 si erg 
3.5 M @ 
0.24 cm -3 
0.003 cm" 3 
8 

20 mG 

■5 X 1U 

3 x 10~ 4 
1 


V s (tsn) 
<j(tsn) 
0"s(tsn) 
M a (t S n) 
-Be (t sn ) 

B d {t sn )(= a Bo) 


3745 yr 
1316 kms" 1 
5.2 
3.1 

25 Mq 

4.6 x 10 50 erg 
104 ^G 


Pc/(poV s 2 ) 


0.145 
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Parameter description: The quantities d and R s denote the assumed dis- 
tance and the radius of the source, respectively, Bd is the internal mag- 
netic field strength, as determined from the thickness of observed X-ray 
filaments cf. eq. 1, and _E sn is the total hydrodynamic explosion energy; 
M e j, M s (t sn ), N S (R S ) and N g (r = 0) are the ejected mass, the swept- 
up mass, the circumstellar gas number density at the SNR shock, and 
the number density at the centre, respectively; k is the power law index 
of the ejecta velocity distribution; Bo is the assumed amplified magnetic 
field strength in the upstream region of the shock precursor, while r\ and 
_R" op denote the assumed proton injection rate and energetic electron-to- 
proton ratio, respectively; t sn is the calculated age of the SNR; V s (t sn ), 
a(t S n), <7 s (t sn ), E c (t sn ), and Bd(i sn ) are the resulting values of the 
subshock velocity, the total compression ratio, the subshock compres- 
sion ratio, the total nonthermal energy, and the downstream magnetic 
field strength, respectively. Finally, P c and po = m p N s (R s ) denote the 
postshock pressure of accelerated particles and postshock mass density 
of the gas, respectively. 



2007a). This confirms the shell-type nature of the source, al- 
though confusion with the emission from the Vela SNR pre- 
vented a detection based on the radio data alone. 

In X-rays, a central diffuse source of ~ 9' x 14' ex- 
tension is seen NW of the geometrical centre of the SNR 
(ISlane et all 1200 U IBecker & Aschenbachl 120021: [Becker et all 
120071) . While ISlane et al.1 (l200ll) argue that the hard spectrum 
of the s ource seen wi t h ASC A hints at a Pulsar Wind Nebula 
(PWN), IBecker et all d2007l) argue, based on XMM-Newton 
data, that the source is soft and reject a PWN nature. Hence, 
there is currently no agreement whether RX J0852. 0-4622 is a 
centre-filled, composite SNR or not. Spatially-integrated fluxes 
would only marginally be affected. There is however the possi- 
bility that such a central PWN would slightly infl uence the VHE 
7-ray radial profile (see Aharonian et all l2007al and Fig.|7j. 
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3.1 .2. Central compact object and pulsar association 



of 



the shell of RX J0852.0- 
C XOU J08520 1 .4-46 1753 
_ 12001b iKargaltsev et al.L 
12002 : Mereghetti, 2001; Mereghe ttietail 120021) , a central com- 



At the geometrical centre 
4622 lies the X-ray po i nt source 
( Aschenbach et all 1 19991: iPavlov etal 



pact object (CCO) similar to that detected in the centre of 
CasA. CXOUJ085201. 4-461753 might be a neutron star, but 
the nature o f the object and the possibly associated compact 
Ha nebul a (Pellizzoni et al., 2002J) is sti l l under debate (e.g . 
iRevnoso et all l2006t iMignani et all 120071: iBeckeretail |2007|) . 



No X-ray pulsations h ave been detected (Karg attsev et all 
2002; lBecker et alll2007l) . The association of CXOU J085201.4- 
461753 with RX J0852.0-4622 is nevertheless suggestive and, if 
true, would allow conclusions on the nature of the progenitor of 
RX J0852. 0-4622. Since the scenario discussed in this paper im- 
plies a core collapse SN, it does not exclude an association of 
CXOU J08520 1 .4-46 1 753 with RX J0852.0-4622. 

In the literature also the possibility is discussed that 
RX J0852.0-4622 i s as sociated with PSR J0855-4644 
(Redman & Meabur ijl 120051) . This appears unlikely how- 
ever, given the implicati ons of this association on di stance and 
age of RX J0852.0-4622 (iRedman & Meaburnll2005l) . 

3.1.3. Non-thermal X-rays 

The soft X-ray emission from RX J0852. 0-4622 is heavily con- 
fused by thermal emission from the Vela SNR. While the emis- 
sion from the Vela SNR seems consistently constrained to two 
thermal components (T12 = 0.05, 1.2 keV) with an absorption 
column density of about 10 20 cm~ 2 dLu & Aschenbachl 2000), 
the surface brightness and temperatures of these components 
are variable enough to prevent a clean subtraction of the Vela 
SNR in the X-ray spectra of RX J0852 . 0-4622 dSlane et all 
I200U llvudin et all 120051: lAharonian etafl l2007al) . In the spec- 
tra above ~ 1 keV, the non-thermal emission from RX J0852.0- 
4622 nevertheless clearly dominates. We adopt the common in- 
terpretation that this component is synchrotron emission from 
relativistic electrons and use two estimates of the total X-ray 
synchrotron flux from RX J0852. 0-4622 in our modelling (see 
Figs. [3] |5j: The first (lower) estimate was derived using the av- 
eraged powe r-law spectra deriv ed from the three brightest shell 
components dSlane et all 12001 1) , which we scaled up to match 
the total flux meas ured with ROSAT in the soft X-ray band 
dAschenbachl [l998). The second (up per) estimate is a reana lysis 
of the ASCA data set as presented in Aharonian et all d2007al) . 

High-resolution imaging of this synch r otron X-ray com- 
ponent with XMM-Newt on (llvudin et all 120051 120071) and 
espec ially with Chandra (Bam ba"et all 120051: iPannuti et all 
|2004|) permits the derivation of synchrotron cooling times and 
therefore of an estimate of the effective magnetic field (e.g . 
Bere zhko et all l2003at iBerezhko & Volkl l2004at iBamba et all 
2005). In order to test our model we will compare the B -field 
derived in this manner to the field value that is required to fit the 
spatially integrated synchrotron data (see section 4). 

3.1.4. Thermal X-rays 

The detection of thermal X-ray emission could help to con- 
strain the gas density and/or its spatial distribution in the SNR. 
However, the interpretation of the X-ray spectra is impeded by 
the strong background emission from the Vela SNR; the latter is 
dominated by low-temperature (< 1 keV) X-ray emission. After 
subtraction of such emission the remaining emission (dominant 



at higher energies) has been attribute d either to gas at highe r 
temperatures from RX J0852.0-4622 (lAschenbach et all 11999). 
or to synchrotron emis sion plus a very small thermal contribu- 
tion dSlane et alll2001l) . The latter conc lusion has basically been 
taken over in the recent H.E.S.S. paper (Aha ronian et al ., 2007a) 
that also contains a re-analysis of the ASCA data, even though 
the X-ray line features exhibited in the spectrum below 2 keV 
might be associated with either RX J0852. 0-4622 or the Vela 
SNR, or to both. Similarly the residuals of the spectra to the 
ASCA data visible around 1 keV are suggestive that another 
component, which might originate from RX J0852. 0-4622, is 
needed. The possible detec tion of a thermal e mission component 
from RX J0852.0-4622 bv lUchivamald2008l) is to be mentioned 
here again. 

These are complex possibilities. The calculation of the ther- 
mal emission for the concrete case of RX J0852. 0-4622 is 
compounded by the wind-bubble plus swept-up-shell geome- 
try. Standard methods for calculating the thermal emission from 
SNRs approximate the configuration either by a plane shock ge- 
ometry or by a classical Sedov solution, and they neglect the ex- 
istence of the accelerated particle component. Neither of these 
approximations is appropriate for the case of RX J0852. 0-4622 
in its present phase. The plane approximation disregards the adi- 
abatic gas cooling in the interior and the overall dynamic evolu- 
tion of the system. The classical Sedov solution implies a uni- 
form circumstellar medium. In the extreme case the radiatively 
cooled wind bubble shell is even denser - and therefore even 
thinner - than assumed in subsection 3.2.2. and it might only 
be reached recently before the present epoch. Then collisional 
electron heating has had little time to operate until now. 

Another important aspect is the modification of the SNR 
shock by the accelerating CRs. As mentioned before, for a 
SNR shock propagating into a uniform medium with a uni- 
form magnetic field this implies that the larger part of the 
shock surface corresponds to a quasi-perpendicular shock with 
a stro ngly reduced injection of nuclear particles dVolk et all 
|2003|F I Suprathermal injection of ions is only possible in the 
quasi-parallel shock regions. If the spatial scales of the quasi- 
perpendicular regions are large enough, then the cross-field dif- 
fusion of the highest-energy particles, accelerated in the mag- 
netic flux tubes delineating the neighboring quasi-parallel shock 
regions, does not reach deeply into these quasi-perpendicular re- 
gions. In the corresponding magnetic flux tubes there are no nu- 
clear particles to be accelerated, there is no magnetic field am- 
plification, and the shock remains unmodified there. This means 
that in the quasi-perpendicular regions the shock dissipation and 
therefore the gas heating occurs in a locally unmodified shock 
with the overall shock speed, leading to a correspondingly high 
gas temperature and high thermal emission. In the case where 
a radiatively cooling shell of a wind bubble is the obstacle for 
the SNR expansion, the situation may be different. MHD insta- 
bilities and the radiative cooling of such a shell probably break 
it into many small regions with strongly varying field direc- 
tions. Then the spatial scales separating the flux tubes originating 
from quasi-perpendicular regions from those of quasi-parallel 
shock regions may become small enough that cross-field diffu- 
sion can also fill the quasi-perpendicular flux tubes with acceler- 
ating particles and then particle acceleration plus magnetic field 
ampl ification occur practical l y eve rywhere over the shock sur- 
face dVolk et all 120081 IVolkL 120081) . In the extreme this implies 
shock modification over the entire shock region and thus a re- 



5 For earlier discussions of this question in a more general context, 
see lEllison etall <1995h : lMalkov & V61kl d!995h . 
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duced gas heating due to the subshock dissipation only. The en- 
hanced overall acceleration efficiency then also demands a lower 
density of the thermal gas for a given hadronic 7-ray flux, and 
thus a lower thermal emission. The low swept-up mass in a low- 
density wind bubble in addition lowers the overall thermal emis- 
sivity compared to that of a classical Sedov remnant with the 
same upstream gas density at the shock at the present epoch. 

In section 4.4 a rough estimate of the resulting thermal X-ray 
emission will be given, based on the emission from a classical 
Sedov solution in a uniform ambient medium. According to this 
estimate the thermal emission of soft X-rays at 1 keV is larger 
than the corresponding nonthermal X-ray emission. However, 
the error in this estimate is not known and is likely to be quite 
large. Therefore an uncertainty remains which we cannot resolve 
at this point. 

3.1 .5. X-ray morphology absorption and relation to CO 
data 

Absorption of soft X-rays by neutral hydrogen can be u sed to 
put constraints on the source distance. ISlane et ail d2001l) used 
CO data of the Vela Molecular Ridge (VMR) to reject a dis- 
tance of RX J0852. 0-4622 of more than 1-2 kpc, based on the 
lack of strong X-ray absorption variation that should have been 
detected across RX J 0852.0-4622. Th i s is br oadly in agreement 
with the inference by Morigu chi et al.l d2001l) that there is an an- 
ticorrelation of X-ray emission and molecular gas traced in CO, 
especially with regard to the VMR at a distance of 1 — 2 kpc. 

In the 2 — 10 keV band, that is expected not to exhibit absorp- 
tion, an X-ray power-law spectrum can be derived. Assuming 
this power-law to continue down to 0.7 keV. [Slane et all f2001 ) 
have then derived an absorption column density of (4.0 ±1.8) ■ 
10 21 cm" 2 for RX J0852.0-4622. Since this column density is 
much larger than the one towards the Vela SNR, they concluded 
that RX J0852. 0-4622 should be at a much larger distance than 
the Vela SNR. 

3.1 .6. Radio spectrum 

We assume that the radio emission is due to synchrotron radi- 
ation. There is no good radio spectrum available fvelaj61.dvior 
the entire remnant. We use the differential flux values given by 
iDuncan & Greenl d2000l) at 2.42 GHz and 1 .40 GHz, the errors of 
which are representing the uncertainty of the background level. 
The spectral index between the two bands has quite a large error 
{a = 0.4 ± 0.5), but for the north-western rim a better spectrum 
{a = 0.40±0.15) could be derived. If this va l ue is re presentative 
for the entire remnant, as IDuncan & Greenl d2000t) argue, then 
this index is somewhat harder than what we expect from a mod- 
ified SNR shock environment, though still compatible within a 
2 er error range. 

3.1 .7. Gamma-ray and X-ray line emission from 
radioactive 44 Ti decay 

The 44 Ti production in a SN explosion depends on progeni- 
tor star mass and explosion-t ype, with a y i eld spa nning two 
orders of magnitude (see, e.g. iRenaud et all [2006a, and refer- 
ences therein). Becaus e of the short lifetime of ~ 80 years (e.g. 
IWietfeldtetafl fl999h . the mere detection of hard X-ray and 
7-ray de-excitation lines at 69.7, 78.4, and 1157keV from the 
44 Ti radioactive decay products can be used to significantly con- 
strain the SNR age. So far, however, these lines have only been 



detected unambiguously from Cassiopeia A (ag e presumably 
~ 33 years), with COMPTEL o nboard CGRO dlvudin et all 
1 1994 . PDS onboard BeppoSAX dVinketall 1200 lb. an d the 
ISGRI imager onboard INTEGRAL dRenaud et alll2006d) . 

For RX J0852. 0-4622, the situation is unfortunately unre- 
solved. From COMPTEL data, the detection of a 7-ray line 
at 1163 ± 16keV, consistent with the 44 Ti 7-ray decay line, 
was cl aimed and predom i nantly attributed t o RX J0852.0- 
4622 dlvudin et all 119981: lAschenbach etafl Il999h . Using 
this 44 Ti line flux, and using a rather high shock velocity, 
lAschenbach et al.l ([1999) derived an age of 680 ± 100 yrs, and a 
distance of 200 pc. Such 7-ray data also suggest a core collapse 
SN event. 

However. ISchonfelder et al.l d2000t) pointed out that the sig- 
nificance of the COMPTEL 44 Ti result is only marginal. And 
the COMPTEL result could so far not be confirmed with the 
INTEGRAL instrume nts SPI and ISGR I A S PI upper limit is 
so far unconstraining (von Kienli n et al ., 2005). Under a point- 
source assumption, the ISGRI non-detection of t he 78keV line 
would be in conflict with the COMPTEL result dRenaud et all 
l2006al) . Ho wever, an extended-so urce analysis has not yet been 
performed dRenaud et all l2006bh . 

Therefore a nearby and rather recent event is not decisively 
excluded from these specific observations, even though the ar- 
guments for it are rather weak. 

3.2. Model parameters for a SN explosion in 1 kpc 
distance 

As discussed in the previous section, the main observational mo- 
tivation to locate RX J0852. 0-4622 at a distance of ~ 1 kpc 
is the much larger column density in neutral hydrogen derived 
from the X-ray spectrum of RX J0852. 0-4622, compared to the 
values for Vela SNR. 

3.2.1 . Reasons for a wind bubble scenario 

The lack (or low level) of thermal X-ray emiss i on an d the strong 
X-ray synchrotron flux already led lSlane et alj d2001l) to the con- 
clusion that RX J 0852.0-4622 coul d be ev olving into a wind 
bubble. Similarly, IDuncan & Greenl ([2000) argue that the un- 
usual radio properties of RX J0852. 0-4622 (bipolar shell mor- 
phology, low surface brightness) could be explained if the SNR 
has so far evolved in a low-density region. 

We can use the non-thermal X-ray plus 7-ray emission from 
RX J0852. 0-4622 to put these arguments for a wind bubble sce- 
nario on a more firm footing, anticipating for the moment that 
the VHE 7-ray emission is indeed dominated by hadronic emis- 
sion. As we shall see in section 4, this latter conclusion follows 
from the fact that with the assumption of a significantly ampli- 
fied magnetic field it is possible to fit all the nonthermal spectra 
as well as the morphology in nonthermal X-rays and gamma- 
rays - and that this amplified field does not exceed that derived 
from the observed X-ray filaments. Such amplification is how- 
ever only possible by the nonthermal pressure of the nuclear par- 
ticles. 

Arguments simil a r to t hose given below have been used 
in iBerezhko & VoUd d2006l) for the comparable case of SNR 
RX J 17 13.7-3946 to which we refer the reader here. 

In order to yield the observed nonthermal X-ray luminosity 
in the case of a uniform IS M the shock speed s hould b e suffi- 
ciently large, V s > 10 3 km/s dBerezhko & Volkl l2004bl) . Given 
that the SNR would already be in the Sedov phase the observed 
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size R a ~ 20 pc, corresponding to a distance of 1 kpc, would 
lead to the age constraint t sn w 0.4i? s /X<3 X 10 3 yr. For a 
typical SN type la explosion energy _E sn = 10 51 erg and ejected 
mass of 1.4M Q this would then imply a very low ISM number 
density Nh 45 10 -2 cm -3 . On the other hand, the peak TeV 7- 
ray luminosity, achieve d during SNR evolution from a type la 
event, roughly scales as dBerezhko & Volklll997h : 



e 7 K y (e 7 ) w 150 



1 cm 



1 kpc 



eV 



(3) 



for / re = 1. Here F 7 (e 7 ) is the integral flux of 7-rays with en- 
ergies greater than e~, This expression shows that for N-r < 
10~ 2 cm -3 we would have to expect an energy flux e 7 i r ' 7 (e 7 ) < 
1.5 eV/(cm 2 s). This is an order of magnitude smaller than the 
observed flux (see section 4). 

Therefore the nonthermal observations make it clear that 
SNR RX J0852. 0-4622 can not correspond to a type la event, 
if the source distance is as large as d — 1 kpc. As a consequence 
we shall consider a core collapse SN event. 

3.2.2. Wind bubble parameters 

The progenitor stars of core collapse SNe that significantly mod- 
ify the density of their environment are massive main-sequence 
stars with initial m asses M; > 15M which have intense winds, 
e.g. (lAbboltH982l) . In the mean, during their evolution in the sur- 
rounding uniform ISM of gas number density po — m p NisM> 
they create a low-density bubble, surrou nded by a shell of 
swept-up and compresse d ISM of radius dWeaver et all 119771; 
IChevalier & Liangfl l989) 



R 



sh 



0.76(0.5MO 3 /p ) 1/5 , 



(4) 



where M is the mass-loss rate of the progenitor star, V w is the 
wind speed, and t w is the duration of the wind phase. The values 
of these parameters are given in tabular form bv lChevalierl(ll982l) 
in terms of Mi. 

In order to determine the SNR shock dynamics inside the 
shell we model the gas number density distribution in the bubble 
and in the shell in the form (e.g. lBerezhko & Volkl, 12006): 



N g = N h + {r/R sh ) 3 ^-^N sh , 



(5) 



where N s h = <j s IiNism is the peak number density per hydrogen 
atom in the shell, N\, is the gas number density inside the bubble, 
typically very small compared with the shell density, and er s h = 
Nsk/NisM is the shell compression ratio. We note that as a result 
of radiative cooling the compression ratio er s h can exceed the 
classical upper limit of 4. The possibility of a very thin wind 
bubble with cr s h ^S> 4 has been indicated in subsection 3.1.4. 
above. 

The mass of the bubble 



M b = (47ri? 3 h /3)mpiV b 



(6) 



is rather small, Mb < M Q , in the case of moderate progenitor 
masses Mi < 20M Q , whereas the shell mass 



Here we use the gas number density distribution N g (r) 
p(r)/m p in the form 



Av 



{0.003 + 0.24[r/(17.5pc)] 12 } cm 



(8) 



M sh = 4nN H hm 



sh m p J drr (r/R sh ) 




3(<T-1) _ 



{4irR 3 sh /3)N ISM m 



is many hundred solar masses. Therefore, during SNR shock 
propagation through the bubble, only a small fraction of its en- 
ergy is given to gas of stellar origin. The main part of the explo- 
sion energy is deposited in the shell. 



which fixes the gas density at the present shock radius and pro- 
vides a consistent fit for all existing data for SNR RX J0852.0- 
4622. Together with Eq.(4), this relation also connects the exter- 
nal density Atsm with the progenitor mass. 

Such a distribution corresponds to a bubble with cr s h = 5 and 
28 < i? s h < 32 pc created by the wind of a main-sequence star 
of initial mass 15Af < Mi < 20M Q in the surrounding ISM of 
hydrogen number density 1 1 < Nj$m < 49 cm -3 , respectively 
(IChevalier & Lianglll989l) . It implies that this bubble is located 
inside a region of dense gas. 



3.2.3. Further parameters that determine the SN 
evolution and CR acceleration 

We shall use the SNR parameters E sn = 1.3 x 10 51 erg, M cj = 
3.5M , and k — 8 which, as will be shown below, give a good 
fit for the observed SNR properties. 

We shall also use an upstream effective magnetic field value 
Bq = 20 pG, which is required to provide the observed syn- 
chrotron flux in the radio and X-ray bands (see below). Such 
a value of Bq is significantly higher than a merely MHD- 
compressed dense-gas magnetic field in the inner part of the 
shell, at densities N K < 0.24 cm~ 3 . 



4. Results and discussion 

In this section we shall discuss the physical characteristics of 
the wind bubble scenario in detail, and compare them with the 
observations. 

The calculated dynamical characteristics of the SNR are 
shown in Fig. 1. From Fig. la one can see that for the assumed 
distance of 1 kpc the calculation fits the observed SNR size 
R s w 17.5 pc at the age t sn — 3745 yr. 

To fit the observed synchrotron and 7-ray spectra (see be- 
low) we assume a proton injection rate 77 = 10~ 3 . This leads to 
a moderate nonlinear modification of the shock which at the cur- 
rent age of t sn — 3745 yrs has a total compression ratio a xs 5.2 
and a subshock compression ratio er s w 3.1 (Fig. lb). All param- 
eters used and the resulting model properties are summarized in 
Table 1. 

For its adopted density, the wind bubble contains only a 
small amount of gas Mb w 0.3M Q . Therefore the SN shock 
deposits only about 20% of the explosion energy during the ini- 
tial 1000 years of propagation through the bubble, as seen in 
Fig. 1. However, up to the current epoch the SN shock has al- 
ready swept up a considerable mass M sw m 25M Q . Therefore 
the ejecta have transformed about 85% of their initial energy into 
gas and CR energy (Fig. lc). The acceleration process is there- 
fore characterised by a high efficiency under the assumption of 
spherical symmetry: at the current time about 35% of the explo- 
sion energy have been transferred to CRs, and the relative CR 
energy content E c continues to increase to a maximum of about 
0M>E sn in the later phase (Fig. lc), when particles start to leave 
the source. 

Therefore, in absolute terms the CRs inside SNR 
RX J0852.0-4622 already contain 



E c « 0.35£ s 



4.6 x 10 50 



erg. 



(9) 
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Fig. 1. Model parameters as a function of time: (a) Shock ra- 
dius R s and overall shock speed V s \ (b) total shock (a) and 
subshock (cr s ) compression ratios; (c) ejecta (E e y) and CR {E c ) 
energies in spherical symmetry. The vertical dotted line marks 
the current epoch of SNR evolution. The external gas density 
iV g = pg/rrip = 1.4-ZVh is also shown in panel (a). 



The volume-integrated (or overall) CR spectrum 

/>oo 

N(p, t) = 16ttV / drr 2 f(r,p,t) 



(10) 



has, for the case of protons, almost a pure power-law form 
N oc over a wide momentum range from 10~ 2 m p c up to 
the cutoff momentum p max ~ 6 x 10 5 ra p c, corresponding to a 
cutoff energy of w 5.6 x 10 14 eV. (see Fig. 2). This value p mSLX 
is limited mainly by geometrical factors, which are the finite size 
and speed of the shock, its decelerati on and the adiab atic cool- 
ing effect in the downstream region dBerezhkol fl996T) . Due to 
the shock modification the power-law index slowly varies from 




l°g(p/m p c) 

Fig. 2. The overall (volume-integrated) CR spectrum as a func- 
tion of particle momentum. Solid and dashed lines correspond 
to protons and electrons, respectively. 
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Fig. 4. Calculated nonthermal Bremsstrahlung (NB, dash-dotted 
line), inverse Compton (IC, dashed line), and 7r°-decay (solid 
line) 7-ray spectral energy distributions as functions of pho- 
ton energy e 7 for the high-injection, high -field model. The 
observed H.E.S S. (lAharonian et all l2007al) and CANGAROO 
(lEnomoto e t al., 2006) 7-ray fluxes are shown as well. 



7 = 2.4 at p ^ m p c to 7 = 1.9 at p ~ 100m p c. The shape 
of the overall electron spectrum N c (p) deviates from that of the 
proton spectrum N(p) at high momenta p > p\ w 350m p c, as 
a result of the synchrotron losses in the downstream region with 
a magnetic field strength £?d ~ 100 fiG which is assumed uni- 
form (i?d = B2 = o~Bq). Therefore within the momentum range 
pi < p < pf nax the electron spectrum is considerably steeper 
N c oc p~ 3 . 

Specifically, the synchrotr on losses become imp ortant for 
electron momenta greater than (Berezh ko et al.Ll2002l) 



Pi 



1.3 



10 s yr\ {10 pG 



f 



(11) 



Substituting the SN age t — 3745 yr into this expression, we 
have p\ w 350m p c, in agreement with the numerical result (see 
Fig. 2). 

The maximum electron momentum can be roughly estimated 
by equating the synchrotron loss time with the acceleration time. 
This gives (e.g. BerezhkoltaLl|2002|) 



Pn 



6.7 x 10 4 



V s 



10 3 km/s 
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Fig. 3. Calculated broadband spectral energy density of RX J0852. 0-4622 , as function of photon energy e 7 . In the 7-ray region 
the solid line shows the ir° — deca y emission and the dash ed line denotes the inverse Compton (IC) emission. Radio fluxes were 
obtained with the Parkes telescope dDuncan & Greenll2000l) . For th e X-ray synchrotro n flux a lower boundary is given by the sum 
of the ASCA fluxes from the brightest parts of the SNR, as given b y lSlane et alj d200ll) . scaled up to match the total remnant's flux 
measured with ROSAT in the soft X-ray range (lAschenbachill998h . The upper boundary comes from a re-analysis of the total ASCA 
data from the remnant, a s given in lAharonian et al.l d2007al) . TeV data are from CANGAROO ( Eno moto et all I2006I) and H.E.S.S. 
( lAharonian et alll2007al) . 
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Fig. 5. Calculated synchrotron emission flux as a function of 
frequency v. The solid line corresponds to the high-injection, 
high-field model, the dashed line corresponds to a hypothetical 
very low injectio n, and therefore unmod ified low-field model. 
Parkes radio data ( Duncan & Green . 2000) and ASCA X-ray data 



dSlane et all I2001 



Aharonian et al 



the caption of Fig. (3}. 



2007a) are shown (see also 



(cr-1) /10/xG 



a(l + <72) I B[ 



At the current epoch V s ~ 1316 km/s, which leads to a maxi- 
mum electron momentum p c max ~ 10 4 m p c, in agreement with 
the numerical results (Fig. [2]). 



As a result of the shock propagating through the wind shell 
(cf. Fig. 1) the SN shock speed decreases rather quickly during 
the period t > 10 3 yr. Therefore, during previous evolutionary 
phases the shock has produced electron spectra with cutoff mo- 
menta p^ ax larger than at the current epoch. Due to this factor 
the spatially integrated electron spectrum has a relatively smooth 
cutoff (see Fig. 2). Together with the synchrotron cooling this 
gives a good fit of the observed X-ray spectrum (see below). 

The present-day parameters B^ = 104 fiG and K cp « 
3 x 10~ 4 lead to good agreement between the calculated and the 
measured spectral energy distribution of the synchrotron emis- 
sion in the radio to X-ray ranges at the present time (Fig. [3). 
The steepening of the electron spectrum at high energies due to 
synchrotron losses and the smooth cutoff of the overall electron 
spectrum together naturally yield a fit to the X-ray data with their 
soft spectrum. Such a smooth spectral behaviour is achieved in 
an assumed upstream field of 20 fiG (which leads to the above 
downstream field B^)- 

Fig. [3] also shows the calculated 7-ray spectral energy 
distributions due to ir° -decay, IC emission, and nonthermal 
Bremsstrahlung, together with the existing experimental data. 

According to the calculation, the hadronic 7-ray production 
exceeds the electron contribution by more than two orders of 
magnitude at all energies. In detail these 7-ray spectra are shown 
in Fig. [4] For energies e 7 = 1 — 100 GeV the 7-ray spectrum is 
el" 2 , hardening from e 7 = 0.1 — 1 TeV, 



(12) c l° se t0 dFj/de 



whereas starting from e 7 « 1 TeV it has a smooth extended cut- 
off despite the comparatively much sharper cutoff of the proton 
energy spectrum, cf. Fig. [2] 

Note that the 7-ray cutoff energy e" lax ~ 0.1p max c is sen- 
sitive to the magnetic field strength B^, since the proton cut- 
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off m omentum has a dependence p max oc R s VsBd (Berezhko, 
1996). It is clearly seen from Fig. 4 that the calculated spectrum 
fits the H.E.S.S. measurements in an acceptable way, at least up 
to ks 5 TeV. However the four highest-energy points tend to lie 
below the theoretical curve. This can be the result of escape of 
the highest-energy protons during the deceleration of the shock 
in the shell. See section |4~T| for further details. 

The hadronic dominance in VHE 7-ray emission which we 
predict here, could be further investigated in the near future by 
the Fermi instrument in the GeV region. Even though at such 
comparatively low 7-ray energies the 7-ray background from 
the diffuse Galactic 7-rays is quite significant, especially for 
such a large, low- su rface brightness source as RX J0852. 0-4622 
dDrurv et all 1 1994b . Fermi should be able to detect the overall 
very high 7-ray flux from RX J0852. 0-4622. It is therefore to be 
expected that the Fermi instrument will confirm our prediction 
that the spatially-integrated 7-ray spectral energy density at 1 
GeV is only a factor w 1.5 lower than at 1 TeV cf. Fig. 4, as a re- 
sult of the nonlinear modification of the acceleration process. If 
the nonlinear modification in the wind bubble is less strong than 
assumed here, then this difference in the spectral energy density 
between 1 GeV and 1 TeV should be even smaller. 

In Fig. 5 we separately present the differential synchrotron 
spectrum, produced at the current epoch. For comparison we 
also show a synchrotron spectrum, which would correspond to 
a hypothetical acceleration scenario in which the proton injec- 
tion rate is taken so small (77 = 10~ 5 ) that the accelerated nu- 
clear CRs do not produce any significant shock modification or 
magnetic field amplification; the value Bq = 5 pG is used in 
this case. There are two small but distinct differences between 
the synchrotron spectra that correspond to these two scenarios. 
The high-injection, high-field scenario leads to a steep differen- 
tial radio frequency spectrum S v oc v~ a with power law in- 
dex a « 0.7, whereas for the unmodified, low-field scenario 
a = 0.5. Unfortunately, the low quality of the existing ra- 
dio data does not allow us to distinguish these two scenarios 
in the radio range, in order to conclude from the radio spec- 
trum alone whether or not we deal in the case of RX J0852.0- 
4622 with efficient CR acceleration leading to a significant shock 
modification and magnetic field amplification. The essentially 
different behaviour of these two spectra at X-ray frequencies 
around v = 10 18 Hz demonstrates on the other hand that in 
the case of strong CR production and amplified magnetic field 
i?d ~ 100 /iG the spectrum S v (v) naturally exhibits a smooth 
cutoff consistent with the experiment. In the simple, unmodified 
low-field case the spectrum S v {v) has too sharp a cutoff to be 
consistent with the experiment. 

It is noted that the X-ray flux represented in Figs. 3 and 5 
comes from two different analyses of the X-ray flux (see also 
section 13. It : the lower boundary was derived by summing up 
the A SCA fluxes from the brightest parts of the SNR dSlane et all 
|2001|) and scaling the re sult up to ma t ch the total SNR's flux as 
measured with ROSAT (lAschenbachl ll998a. The upper bound- 
ary - a factor of two higher - comes from a r e-analysis of the 
total A SCA data from the remnant, as given in lAharonian et alj 
(I2007al) . 

The properties of small-scale structures of RX J0852. 0-4622 
seen in X-rays furnish even stronger evidence that the magnetic 
field inside the SNR is indeed considerably amplified. As in the 
case of other young SNRs (e.g. SN 1006, Cassiopeia A, Tycho's 
SNR) Chandra shows very fine filamentary structures in non- 
thermal X-rays in the very outer part of the remnant. The thinnest 
filament detected bv lBamba etafl d2005l) has an angular thick- 
ness Ai/> « 38" in the radial profile of the X-ray emission in 
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Fig. 6. The projected radial profile of the X-ray synchrotron 
emission for e„ = 1 keV. The abscissa is normalized to the 
shock radius and corresponds therefore to an angular scale. The 
solid line corresponds to the high-injection, high-field model, 
the dash-dotted line corresponds to a hypothetical low-injection, 
low-field model. 

the 2-10 keV range. In order to find out whether this type of 
structure is consistent with our theory we present in Fig. 6 the 
projected radial profile 

J(e,p)oc J dxq(e, r = \/ p 2 + x 2 , x), (13) 

calculated for the X-ray energy e„ = 1 keV. The abscissa in 
Fig-H] (and correspondingly also in Fig.|7]i is scaled to the ra- 
dius of the remnant, i.e. to 17.5 pc, as derived from the angular 
radius of the remnant of ~ 1°. Here q(e,r) is the luminosity 
in the nonthermal emission with photon energy e. The integra- 
tion is performed along the line of sight. One can see that the 
theory predicts the peak of the emission just behind the shock 
front with a thickness Ap/R s « 10~ 2 which corresponds to an 
angular width ss 36". It corresponds reasonably well to the 
Chandra observation. At the same time Fig. 6 also shows that the 
unmodified low-field solution cannot explain the filament struc- 
ture seen with Chandra. 

It has been already demonstrated for other young SNRs 
dBerezhko et al.Ll2003at berezhko & Volkll2004al) that the mea- 
sured width of the projected radial profile of the nonthermal X- 
ray emission gives the possibility to determine the internal mag- 
netic field strength according to Eq.(l). 

Substituting into this equation I2 = L/7 = 8.2 x 10 16 cm, 
a = 5.2, V s = 1316 km/s, and v = 3 x 10 17 Hz (i.e. X-ray en- 
ergy e„ ss 1 keV), we obtain w 139 /iG, which agrees within 
30% with the value B^ = 104 /iG used in our spectrum calcu- 
lation. Such a difference in B& corresponds to the uncertainty in 
the field determination in the other objects analysed up to now. 

The magnetic field amplification is driven by the gradient of 
the (nuclear) CR pressure upstream of the outer shock, and we 
can check whether RX J0852. 0-4622 belongs to the class of ob- 
jects that fulfil Eq. 2. In the present case of d = 1 kpc we have 
P c « 0.15poV 2 , where po = N g (R s )m p is the ambient gas den- 
sity at the current shock front position. Substituting N g (R s ) — 
0.24 cm" 3 and B = 20 pG we have B 2 /{8irP c ) w 6.5 x 10" 3 , 
in rather good agreement with the average number in Eq.(2). 
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Fig. 7. The 7-ray emissivities for the 7-ray energy e 7 = 1 TeV 
as a function of projected, normalized radial distance p/R s . The 
calculated radial profile is represented by the solid line. Data 
points are from the Northern part of RX J08 52. 0-4622 as mea- 
sured with H.E.S.S. (lAha ronian et al., 2007a), with an analysis 
point spread function of Gaussian width 0.06°. The dashed line 
represents the calculated profile convolved with the same point 
spread function. 



The line-of-sight integrated 7-ray emission profile as a func- 
tion of projected radius p is calculated for e 7 = 1 TeV and is 
presented in Fig. 7. Due to the large radial gradient of the gas 
and CR distributions inside the SNR, the theoretically predicted 
three-dimensional radial emissivity profile of TeV-emission is 
concentrated within a thin shell of width Ar ~ 0.01i? s . As a 
result of the projection effect, the two-dimensional width is by 
a factor of seven larger than the width of the three-dimensional 
profile, i.e. Ap w 0.1i? s (solid lines in Fig.[7j>. Since the H.E.S.S. 
instrument in addition has a finite angular resolution we present 
as dashed lines in Fig. 7 also the modified radial profile con- 
volved with a Gaussian point spread function with a a p = Ap = 
0.06i? s , corresponding to an angular resolution of <j^, — 0.06°. 

Fig. 7 shows that the expected radial profile of the TeV- 
emission, after taking into account the instrumental angular res- 
olution, is much broader than the intrinsic projected profile 
and is characterised by a minimum-to-maximum intensity ratio 
J™ in /J™ ax w 0.35. The radial profile measured by H.E.S.S. 
compares reasonably well with the theoretical prediction if we 
ignore the two data points in the central region. This might be 
justified because a central compone nt (e.g. a PWN) cannot be 
fully excluded, as was also argued in Ahar onian et alj d2007al) . 

To illustrate the expected time evolution of the nonthermal 
emission we give in Fig. 8 the fluxes of radio emission at the fre- 
quency v = 1 GHz and X-ray emission at the energy e„ = 1 keV, 
as well as the TeV flux, all as a function of time. Since the CR en- 
ergy content still increases, the radio emission increases during 
the next centuries. However, due to the substantial shock decel- 
eration the nonthermal X-ray emission is expected to decrease 
with time at a rate of about 0.04%/yr. The TeV 7-ray emission, 
on the other hand, will increase at a rate of 0. 14 %/yr due to the 
continuing increase of the ambient gas density. Nevertheless, a 
change of the order of one percent in 7 years is probably very 
difficult to measure. Therefore we can not put much practical 
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Fig. 8. The time dependence of the fluxes of the synchrotron 
emission at frequency v = 1 GHz (dash-dotted line), syn- 
chrotron X-ray emission with energy e„ = 1 keV (dashed line) 
and TeV-energy 7-ray (solid line). The fluxes are normalised to 
their values Fo at the current epoch. 

weight on these expectations regarding the secular variation of 
nonthermal fluxes. 

4.1. Escape 

The decrease of the shock speed during SNR evolution di- 
minishes the maximum energy to which particles can be accel- 
erated during any given phase. This has the well-known con- 
sequence that particles already accelerated to a higher maxi- 
mum energy at earlier times can now begin to diffusively es- 
cape from the interior of the remnant, even without any change 
of the scattering mean free path, as it was predicted analyti- 
cally ( Berezhko & Krymskvl Il988l) and confirmed numerically 
dBerezhkoet alj,|j996|). 

Compared with this escape - which is relatively slow in 
our model because the diffusion of CR particles with momenta 
P > Pm(t) higher than the current cutoff momentum p m (t) is 
still described by Bohm diffusion - the actual CR escape is 
presumably much faster. Since the CRs with p > p m (t) de- 
velop a much smaller spatial gradient than those with momenta 
P < Pm{t) (which are still efficiently accelerated by the shock), 
their production of chaotic magnetic fields whose scales are ap- 
proximately equal to their gyroradii goes down quickly with 
time. As a consequence their diffusion coefficient increases con- 
siderably so that they leave the SNR much more rapidly than 
through Bohm diffusion, and then these CRs contribute much less 
to the 7-ray flux of the source than calculated above. The 7-ray 
spectrum corresponding to the particles that remain confined in 
the SNR should therefore drop off faster with energy than given 
by the solid line in Fig. 4. However, the H.E.S.S. data at the high- 
est 7-ray energies are not suggesting that this effect is very sig- 
nificant in the present case. 

Secondly, the field amplification is probably stronger dur- 
ing earlier phases of the remnant evolution , when poV^ is larger 
than at later phases dBell & Luceki |2001|) . Our calculation has 
assumed the effective field strength to be constant. Correcting 
for this, the overall nuclear particle spectrum released into the 
Interstellar Medium by SNRs presumably forms the Galactic 
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CR s pectrum up to an energy of 10 17 eV (Berezh ko & Volkl 
|2007|) . That protons alone should be able to reach energies 
> 10 15 e V was shown for t he case of a representative type 
la SN by iBerezhko & Volkl d2004bl) . This is likely to be true 
also for RX J0852. 0-4622 and, by implication, also for its twin 
RX J 17 13.7-3946, even though their environmental conditions 
are drastically different from those of a type la SN. 

We add here that the above maximum energies of energetic 
nuclei are calculated under the assumption of Bohm diffusion. 
To this extent they are upper limits, given the strength of the 
corresponding amplified magnetic field. 

Except in special cases one should therefore ex- 
pect that escape sets in when th e source becomes older. 
iPtuskin & Zirakashvilil (120031 120051) have in addition argued 
that at late times the damping of the scattering magnetic 
fluctuations should increase the scattering mean free path. The 
decrease of the effective field strength on the one hand, and of 
the scattering strength (through the slowing-down of the shock 
and wave damping) on the other, go in the same direction to 
lower the cutoff energy of the accelerating population. The 
relative importance of these effects probably depends on the 
details of the source, in particular its evolutionary phase. 

4.2. A rough estimate for the thermal X-ray emission 

Using the results for the dynamical evolution of the system one 
can also attempt to estimate the thermal X-ray emission. As ar- 
gued in section 3. 1 .4 this is only possible in an approximate way, 
even if the ejecta emission is disregarded with the argument that 
the ejected mass in the present phase is small compared to the 
swept-up mass. In fact, looking at Fig. la it can be seen that the 
SNR is well past the sweep-up phase and has entered a quasi- 
Sedov phase in the stellar wind shell, i.e. a roughly self-similar 
evolutionary phase, modified by strong particle acceleration rel- 
ative to a purely gas dynamic evolution. 

The approximation used is the following. The bubble case 
is compared with a SNR in a uniform medium in the classical 
Sedov phase without any CR acceleration, making four assump- 
tions (i) the total hydrodynamic explosion energy is the same in 
both cases (ii) the downstream unmodified temperature T s , de- 
termining the overall shock velocity and the thermal emission, 
is the same, (iii) the present gas density upstream of the shock 
is the same, and (iv) the t wo objects are at the s ame distance of 
1 kpc. Then the results of lHamilton et alj d 1983b for the thermal 
X-ray flux from a classical Sedov SNR are used, employing the 
emission measure of the bubble remnant instead of that of the 
classical Sedov remnant with the same five parameters above. 
This means that the X-ray emissivity of the remnant is reduced 
by the ratio i? om = EM\> /EM$ of the emission measure EM^ 
for the bubble solution to the emission measure for the classical 
Sedov solution EM$ which corresponds to a uniform ambient 
gas density N g (R s ). 

The thermally relevant gas temperature in the presence of 
particle acceleration is the temperature T su b downstream of the 
gas subshock. It corresponds to the temperature downstream of 
a shock that is nonlinearly modified by the internal energy and 
the pressure of the accelerated particles. 

The CR-modified shock, i.e. the precursor-subshock system, 
is approximated by a plane parallel structure, which implies that 
the precursor size is small compared to the SNR shock radius. 
This is generally the case. Then conservation of mass, momen- 
tum and energy fluxes permits in a straightforward way to cal- 
culate the ratio between the actual downstream gas temperature 



T su b and the downstream gas temperature T s that would obtain 
without particle acceleration: 

Tsub [(a s -l)(7 + l) + 2]( 7 +l) 2 
T s 4 7 ( 7 -l) ( r 2 ■ ^ } 

Here a and a s are the total and the subshock compression 
ratio, respectively, as given in Fig. 1 . The specific heat ratio 7 for 
the nonrelativistic thermal gas may be taken as 7 = 5/3. For the 
present phase of RX J0852. 0-4622 , i.e. a w 5.2 and a s rj 3.1, 
one t hen obtains T su b w 0.45r s . 

lHamilton et alJdl983h give the gas temperature T s by 

T s = 10 7 (U s /839 km s -1 ) 2 K (15) 

in terms of the overall shock velocity V s , assuming a highly ion- 
ized system. For RX J0852. 0-4622 the latter quantity is given 
in our Fig. la. Using the present value V s = 1316 km/sec, the 
foregoing equations result in T s w 2.5 x 10 7 K and T su b ~ 
1.1 x 10 7 K, corresponding to a thermal energy of about 1 keV. 

For the emission measure EM we have EM = 
drr 2 N 2 {r), if we disregard irrelevant numerical factors. 
In order to calculate R cm for the CR-modified shock the over- 
all shock is approximated by a shock in a thermal gas with 
an adiabatic index 7 cr such that (7 cr + l)/(7 cr — 1) = <r, 
where a denotes, as before, the total compression ratio of the 
CR-modified shock. The corresponding solution for the g as dy- 
namic quantities is approximated as being self-similar cf. lSedovl 
(1959). If the ambient density profile has a radial dependence 
in the form of a power law N g (r) = No(r / i?o)' 3 , then the 
self-similar density distribution in the SNR, downstream of the 
shock, can be further approximated by a power-law profile, with 
the same swept-up mass as the full solution. It has then the 
form iV g (r) = aN {r / R^f 3 ' with j3' = 3(cr - 1) + /3a, where 
No = N g (R s ). Such a density profile gives the emission mea- 
sure EM = N^Rla 2 /{2(3' + 3). 

The classical gas dynamic Sedov solution for a SN explosion 
into a uniform medium has (3 = and er = 4. Using the same 
values for No for both density profiles leads to R s 12.51 pc for 
the Sedov case and to i? cm w 0.63 for the bubble case (3 = 12 
and g = 5.2. The swept-up mass in the Sedov case is Ms ~ 
50.3 M Q . 

Using the differential therma l X-ray model s pectra dF/ de « 
10~ 4 photons/(keV cm 2 s) from lHamilton et all dl983) (see their 
Fig. 2) for their lgT s = 7.75 and lgT s = 7.25, 77 = N^E sn = 
10 49 erg cm -6 , scaling it according to dF/ de oc i] with a fac- 
tor of sw 3.8 for 77 sw 3.8 x 10 49 erg cm -6 , multiplying it by 
the factor 6 2 [E sn /(10 51 eTg)]- 1 / 2 « 0.65 x 10 4 , as required, 
where 9 ks 86 arcmin for the angular size of the classical Sedov 
remn ant corresponding to the above parameters (lHamilton et al.L 
1983), and multiplying it finally also by the factor R cm , results 
in a thermal spectral energy density e 2 dF/de « 1560 eV cm -2 
sec -1 for e = 1 keV. 

This must be compared with the observed nonthermal X-ray 
energyfluxat 1 keV (see Fig. 3) e 2 dF/de 100eVcm~ 2 sec -1 . 
Therefore, at 1 keV, the thermal flux, calculated in this form, 
comes out to be larger than the nonthermal flux by a factor of 
about 16. 

We note however that the above X-ray energy flux has to 
be considered as a rough upper limit estimate because of the 
following reasons: 

First of all, the SN shock into the bubble and shell interacts 
with a strongly rising gas density profile. Yet for the estimate the 
peak value of the gas density was used. 
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Moreover,in the case of a modified shock the actual gas tem- 
perature T su b is at least by a factor of 2 lower than the temper- 
ature T s used in our estimate. From a general point of view the 
actual thermal X-ray emission is expected to be lower due to 
lower temperature. A rough correction can be done in the fol- 
lowing way. The shock radius corresponding to the Sedov solu- 
tion depends on the relevent parameters according to the relation 
R s oc (EsN/iTsNft)) 1 / 3 . According to this relation the shock 
with the same radius but with two times lower postshock temper- 
ature corresponds to a two times lower explosion energy. Such a 
shock is characterised by the parameter value 77 rj 1.9 x 10 49 erg 
cm -6 that gives a 1 keV thermal X-ray flux which is lower by a 
factor of two compared with the above consideration. 

Third, as it was already mentioned, the calculation in section 
4 was performed for a rather moderate value of the injection rate, 
which leads to a moderate shock modification. Higher injection 
rates can not be excluded given the present data which can lead 
to a considerably higher shock modification with significantly 
lower gas energy and gas temperature. For a comparative ly ex- 
treme position, see the recent paper bv lDrurv et al.l d2008). 

Fourth, for a fixed 7-ray flux F 1 oc NhE$n /d 2 the re- 
quired gas number density scales as Njj oc d 2 with source 
distance d. Consequently a decrease of d requires a decrease 
of the gas density. The thermal X-ray luminosity is quite sen- 
sitive to the gas density and to the total gas energy because 
roughly dF/de oc NfjEsN- A somewhat smaller distance than 
the d = 1 kpc, assumed here, can not be excluded. 

Finally, the bubble solution developed in this paper has quite 
a different temperature profile than any "equivalent" classical 
Sedov solution used to estimate the thermal emission. This arises 
from the approximate uniformity of the total pressure in the SNR 
interior which roughly implies T su b oc N^ 1 . The density pro- 
files are indeed quite different in both cases. 

As a result, we believe that the uncertainties regarding the 
thermal X-ray emission are quite large. This situation also im- 
plies some systematic error in the overall model presented, al- 
though this can hardly change its key properties. In our opinion 
the uncertainty in the thermal emissi on has to be re s olved by 
future work that extends the efforts of Ha milton et al.l d 1983b to 
more general circumstellar density profiles and ejected masses 
for core-collapse supernovae. 

5. Conclusions 

We have argued in the last section that a solution behind the 
Vela SNR, involving the core collapse of a massive star in its 
own wind bubble, describes the available data reasonably well. 
The physical characteristics of this solution are very similar to 
those of the well-known SNR RX J1713.7-3946. In this sense 
RX J0852.0-4622 and SNR RX J1713.7-3946 are twins. 

The calculated 7-ray spectrum has a smooth cutoff at higher 
energies. The presently observable lower limit for the magnetic 
field amplification - in terms of thin shock filaments in hard X- 
rays - is consistent with the one deduced from a theoretical fit 
to the observed synchrotron spectra. If anything, the spectrally 
deduced field amplification is lower than the filament-deduced 
field amplification. Overall, the field amplification is also con- 
sistent with the semi-empirical relation Eq.(2) between the CR 
pressure P c and the magnetic pressure, where P c stems from the 
proton injection rate required to fit the observed 7-ray emission. 
We conclude from these consistency arguments that the obser- 
vational lack of a detailed radio synchrotron spectrum does not 
preclude the determination of a consistent amplified field. 



The magnetic field amplification results in a significant de- 
pression of the density of ultra-high energy electrons and reduces 
the IC and NB contributions of these electrons to less than one 
percent of the 7r°-decay 7-ray emission. 

A remaining uncertainty is connected with the thermal emis- 
sion properties. For the wind bubble configuration this may in- 
deed not be too critical, since gas heating should occur primarily 
at the subshock alone. In addition, the amount of swept-up mass 
is rather small compared to a classical Sedov remnant in a uni- 
form circumstellar medium of equal present preshock density. 

Apart from this uncertainty the main result is the hadronic 
dominance in the 7-ray emission spectrum. As a consequence of 
the nonlinear modification of the shock the spatially integrated 
7-ray spectral energy distribution at 1 GeV is predicted to be at 
best a factor 1.5 lower than at 1 TeV. 
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